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§ : ABSTRACT 

^ . In the frame of the Lindblad theory of open quantum systems, the system of three 

! uncoupled harmonic oscillators with opening operators linear in the coordinates and mo- 

^ I menta of the considered system is analyzed. The damping of the angular momentum and 

^ ' of its projection is obtained. 

^1 1. Introduction 

^ i In the last two decades, more and more interest arose about the problem of dissipation 

O ' in quantum mechanics, i.e. the consistent description of open quantum systems [1-6]. 

O ■ Because dissipative processes imply irreversibility, and, therefore, a preferred direction 

^ I in time, it is generally thought that quantum dynamical semigroups are the basic tools to 

■ introduce dissipation in quantum mechanics. The most general form of the generators of 
^ ■ such semigroups (under some topological conditions) was given by Lindblad [7-9], whose 

. formalism has recently been applied to various physical phenomena, for instance, to the 
^ ! damping of collective modes in deep inelastic collisions [10-15]. An important feature of 
^ I these reactions is the dissipation of energy and angular momentum out of the collective 
degrees of freedom into the intrinsic or single-particle degrees of freedom. 

. !Zi i The damping of energy and angular momentum can be described at different levels of 

^ ' approximations: one may use quantum or classical mechanical methods and microscopical 

■ or phenomenological models. 

In the previous papers [13-15] simple phenomenological models for the damping of the 
angular momentum were studied in the framework of the Lindblad quantum mechanical 
theory. In the present paper, for a special choice of the opening coefficients we obtain ex- 
plicit expressions for the damping of the angular momentum and projection of the angular 
momentum of three harmonic oscillators. As opening operators we use the generators of 
the Heisenberg group, the coordinates qk and momenta pk{k = 1,2,3). In order to keep 
the symmetry of the system we assume that the opening coefficients are symmetric under 
the permutation of the axes. In this way we include the influence of the coupling, due to 
the environment, of the three harmonic oscillators, initially uncoupled. 

The paper is organized as follows: In Sec. 2 we ffist present the Lindblad equation of 
motion in the Heisenberg picture. In Sec. 3 we obtain the damping of the projection of 
the angular momentum for a system of three uncoupled harmonic oscillators. In Sec. 4 we 
analyze the behaviour of the squared angular momentum for the same system. 
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2. The Lindblad formalism in the Heisenberg picture 

In Lindblad's formalism, the usual von Neumann-Liouvillc equations ruling the time 
evolution of closed quantum systems are replaced by the following ones [7-9]: 



dt 
dt 



= L[MA)], (2.1) 



in the Schrodinger and Heisenberg picture, respectively. Here, is the dynamical 

semigroup describing the irreversible time evolution of the open system in the Schrodinger 
(Heisenberg) representation; p is the density operator and A any operator acting on the 
Hilbert space I1{A e B(H)); finally, L{L) is the infinitesimal generator of the dynamical 
semigroup ^t{^t)- 

By using the Lindblad theorem [7-9] which gives the most general form taken by the 
generator L, the Markovian master equation (2.1) takes the form: 

= L{MA)) = ^[hMA)] + ^J2(^;[MA),Vj\ + [V+MA)]V,). (2.2) 

j 

Here, H is the Hamiltonian of the system and the operators Vj , V^^ are taken as polynomials 
of only first degree in the observables qi,q2,Q3,Pi,P2,P3, which represent the hermitian 
generators of the Heisenberg group. Then in the linear space spanned by qk,Pk{k = 1, 2, 3), 
there exist only six linearly independent operators Vj=i^2,...,6' 

3 3 

Vj Cl'jkPk + ^ bjkQk, 

k=i fc=i 

where a^/j, bjk e C withi = 1,2, ...,6. 
Then it yields 

3 3 

k=l k=l 

where a*^, 6*^ are the complex conjugates of ajfc, bjk- 

The coordinates qk and the momenta obey the usual commutation relations (A;, / = 
1,2,3): 

[qk-,pi\ = ifiSki, [qk,qi] = \pk,pi] = o. 

Inserting the operators Vj and V^^ into (2.2) we obtain: 

HA) = ^[H,A] + ^Y^hId^.h - ^akl)iPk[A,Pl] - [A,pk]pi)+ 

k.l 
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2 2 
+i-^DpkPi -0kl){qk[A qi] - [A, qk]qi) - j^D^^p^ {pk[A, qi] - [A,pk]qi+ qi[A,Pk] - [A, qi]pk)- 

-i>^kl{Pk[A, qi] - [A,pk\qi - qi[A,Pk] + [A,qi]pk)}. 
Here we used the following abbreviations: 

DqkQi = Dqiq^ = ^Re{a.l ■ a/), 
DpkPi = ^PiPk = 2^^^^k • b/)) 

h 

Dqupi = Dpidk = ■ hi), (2.3) 

aki = -mk = -Irn{al ■ ai), 
Pkl = -Plk = -Irn{hl - hi), 
Xki = -Im(a.l - hi). 

The scalar products are formed with the vectors a^, hk and their complex conjugates 
a^,b^. The vectors have the components 

a/c = (oifc, ci2k, oefc), 

3. The projection of the angular momentum 

The general Hamiltonian of three uncoupled oscillators is 

H = Ti^^Pl + "^qD- 
^^2mk 2 ^''^ 

The time-dependent expectation values of self-adjoint operators A and B can be writ- 
ten with the density operator p, describing the initial state of the quantum system, as 
follows: ^ 

mA{t) = Ti{p^t{A)), 

and, respectively, 

(JABit) = ^Ti{p^t{AB + BA)). 

In the following we denote the vector with the six components mq^{t),mp^{t),^ = 
1, 2, 3, by m(t) and the following 6x6 matrix by a{t){i,j = 1, 2, 3): 



^PiQj ^PiPj 
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Then via direct calculation of L(qk) and L(p)~) we obtain 

dm 



dt 



= Yin, 



(3.1) 



where 



Y 



From (3.1) it follows that 





-All 


— A12 


— Ai3 


1 

mi 


-OC12 


-ai3\ 




— A21 


— A22 


— A23 


-OL21 


1 

m2 


-0:23 




— A31 


— A32 


— A33 


-Q!31 


-0:32 


1 

mz 




—miuj^ 


P12 


Pl3 


-An 


— A21 


— A31 






—m20J2 




— A12 


— A22 


— A32 


v 


P31 


PS2 




— Ai3 


— A23 


— A33 / 



(3.2) 



(3.3) 



m(t) = M(t)m(0) = exp(ty)m(0), 
where m(0) is given by the initial conditions. The matrix M{t) has to fulfil the condition 

lim M{t) = 0. (3.4) 

In order that this limit exists, Y must have only eigenvalues with negative real parts. 

By direct calculation of L{q}.qi), L{p}.pi) and L{q}.pi + piqk), {k, I = 1,2, 3), we obtain 



da 
'dt 



Ya + aY^ + 2D, 
where D is the matrix of the diffusion coefficients = 1, 2, 3) 



(3.5) 



and Y'^ the transposed matrix of Y. The time-dependent solution of (3.5) can be written 
as 

&{t) = M{t){a{0)-±)M^{t) + ±, (3.6) 

where M{t) is defined in (3.3). The matrix E is time independent and solves the static 
problem (3.5) {da/dt = 0): 

Y± + ±Y^ + 2D = 0. (3.7) 
Now we assume that the following limit exists for t — > 00: 



a{oo) = lim a{t). 

t—^00 



In that case it follows from (3.6) and (3.4): 

(t(oo) = S. 
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(3.8) 



(3.9) 



Inserting (3.9) into (3.6) we obtain the basic equations for our purposes: 



where 



a{t) = M{t){a{0) - d{oo))M^{t) + (7(00), 
f<j(oo) + CT(oo)f^ = -2D. 



(3.10) 
(3.11) 



For the calculation of the matrix M{t) we must diagonalize the matrix Y by solving 
the corresponding secular equation, i.e. det(y — zl) = 0, where z is the eigenvalue and 
/ is the unit matrix. According to (3.2) one obtains an equation of sixth order for the 
eigenvalues z, which can be simply solved only for special examples. In the particular case 
with aki = 0, Pkl = 0, Xki = 0{k 7^ I), the secular equation is obtained as 

[{z + Xiif + ujf] [{z + A22)' + a;i] [{z + A33)' + col] = 0- 

The roots of this equation are 

zi,4 = -All ± Z2,5 = -A22 ± 1(^2, Z3,e = -A33 ± ic^s- (3.12) 

Only positive values of An, A22, A33 fulfil (3.4). Applying the eigenvalues Zi of Y we can 
write the time-dependent matrix M{t) as follows: 

Mmn{t) = J2^rnieMzit)Nr\ 
i 

where the matrix N represents the eigenvectors of Y: 

^ ^ Ymn-^ni — ZiNrni- 
n 

With the relations Mmn{t = 0) = dmn and dMmn{t) / dt\t=o = ^mn and using (3. 3), (3. 10) 
we conclude that the expectation values of the coordinates and momenta decay with 
the exponential factors exp(— Aiit),exp(— A22^) and exp(— Asst) and the matrix elements 
Cmn with the combined factors exp(— 2Aiit), exp(— A22i),exp(— A33t),exp[— (An + X22)t], 
exp[-(Aii + A33)t] and exp[-(A22 + Ass)^]. 

We present here the matrix M{t) only for our special and simple case that the oscil- 
lators are uncoupled. With the roots given in (3.12) we obtain 



M{t) 



( Ml 








Mil 













M2 








M12 













M3 








Mi3 




M21 








Ml 













M22 








M2 







V 





M23 








M3 


) 



(3.13) 



where (/c = 1,2,3) 



Mfe = exp(-Afefet) cosoifet. 



5 



Mik = — - — exp(-Afcfct) sinujkt, (3-14) 
mkUJk 

M2k = -mfca;fcexp(-Afcfct) sinukt. 

This matrix can be used to evaluate d-{t) defined by (3.10). Since we are interested to 
obtain the expectation value of the projection of the angular momentum, we write down 
the following expressions for aq^p^ and (Tq^p^ with M{t) of (3. 13), (3. 14): 

(^qiP2{i) = exp[-(All + A22)t]((o-gip2(0) - CJq,p2(00)) COS Wit COS ^2^+ 

1 



miu}\ 

1712002 



('^PiP2(0) ~ crpiP2(oo)) smLJitcosLJ2t - m2a;2 (0-3152(0) ~ "^9192(00)) cosuJitsmuj2t- 
(crg2pi(0) - crq2pi(oo))sina;itsina;2t) + aq-,p^(oo), (3.15) 

TTllLOi 

(^q2Piit) = exp[-(Aii + A22)t](((7-g2pi(0) - ctqjpi (00)) COS a;it COS 0)2^- 
-mia;i((Tqiq2(0) ~ ^qiq2{'^)) sina;i/; COS a;2t H (c'pip2(0) ~ '^pip2(^)) cosa;i/;sina;2t— 

1712^02 

"^^^^■(^giP2(0) - (Tqip2{oo)) smujitsinuj2t) + aq^p^{oo). (3.16) 



777,20)2 



Similar expressions are found for the other matrix elements of &{t). The matrix elements 
of a{oo) depend on Y and D and must be evaluated with (3.11) or by the relation [11]: 



POO 

a{oo) = 2 / M{t')DM^{t')dt'. 
Jo 



We obtain that 

(^qiP2(oo) = ^g2Pi(oo). (3.17) 

The expectation value of the projection of the angular momentum can be written from 
(3.15) and (3.16): 

<L3it)>=aq,p,it)-aq,p,it). (3.18) 
If the three uncoupled oscillators have the same mass and frequency, then we obtain 

< Ls{t) >= iaq,p,{0) - (7g2Pi(0))exp[-(Aii + A22)t]. (3.19) 

If we consider the open system to be symmetric, then from (2.3)we have An = A22 = A 
and (3.19) becomes: 

< Lsit) >=< Ls{0) > exp(-2At), 

with < Ls{t) > \t^oo ^ if A > 0. This evolution law for the projection of the angular 
momentum is identical to that obtained in [14] by another method. We should like to 
mention that the same evolution law for the projection < Ls{t) > of the angular momentum 
can also be obtained for only two uncoupled harmonic oscillators. 
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4. The squared angular momentum 

In this Section we should hke to discuss the behaviour of the squared angular momen- 
tum < L^(t) > for a system of three independent damped harmonic oscillators. In this 
case the expectation value of the angular momentum can be obtained directly from the 
following expression: 

< L\t) >=< pl{t) >< ql{t) > + < pl{t) >< ql{t) > + < pl{t) >< ql{t) > + 

+ < pl{t) >< ql{t) > + < pl{t) >< ql{t) > + < pl{t) >< ql{t) > - 

< ipiQl+Qm)it) >< iP2q2 + q2P2){t) > < (P292 + 92P2)(i) >< (p393 + Q'3P3)(i) > - 

1 3 
-2 < (PsQs + Q3P3){t) >< {Piqi + qiPi){t) > (4.1) 

The expectation values in (4.1) can be obtained by using the expressions for the centroids 
and variances for the one-dimensional harmonic oscillator obtained in [10]. We have {k — 
1,2,3): 

< Pfe(^) >= (^(pp)fc(*) + f^pfc(^), 

<gfe(0>=^(,,).(0+<(0, (4-2) 

\ < {pkqk+qkPk){t) >= (^{pq)M +(^pAt)(^qkit)- 
For each harmonic oscillator the Hamiltonian is chosen of the form (k=l,2,3) 

= 2m^^'^ ~~2~^^ Y^-^''^'' ^''^''^ 
and the operators in the Lindblad generator can be written in the form 

V> = a\ 'Pk + bl 'qk,i^ 1,2, 

with a\'^\ bf^^ complex numbers. 

In the following we consider, for simplicity, the system of three harmonic oscilla- 
tors with the same mass m, the same frequency u) and the same opening coefficients 
Dqq, Dpp, Dpq and A. In addition, we take ^ = 0. Then, by introducing the expressions 
(4.2) taken from [10] into (4.1), we obtain, after long, but straightforward calculations, the 
following form for the expectation value of the squared angular momentum: 

< L\t) >= pe-^^'+ 

_l_g-2At|^^ sin^ ut + /32COS^ ut + /33siB.Lvt + /34^cosujt + (3^ sina;t cosa;t}-|- 

+l2(oo), 
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where Pj{j = 1, ...,5) are constants. The asymptotic value (if A > 0) is given by 



< (oo) >= (oo)(J(gq), (oo) + (oo)(j(pp), (oo) - 2(j(pq)j (oo)(j(pq), (oo) + {cyclic) 




where [10] 



1 



{{mujf{2X^ + uj^)Dqq + oj'^Dpp + 2muj'^XDpq), 



2(ma;)2A(A2+a;2) 




{{mLofuj^Dqq + (2A^ + uj^)Dpp - 2muj^XDpq), 



2mA(A2 + a;2) 



{-X{mu;fDqq + XDpp + 2mX^Dpq). 



If A > 0, this evolution law shows an exponential damping of the angular momentum of 
the system of three harmonic oscillators, like in [14] . 
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